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Key Points: 


e Correlation between cloud cover and aerosol optical depth is positive at most locations 


except deserts, for all seasons and aerosol types. 


e Increased cloudiness is associated with populations of either smaller or larger aerosol 


particles. 


e¢ Quality assessment flags based on local variability help identifying aerosol populations 


affected by surrounding clouds. 
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Abstract 


This study seeks to help better understand aerosol-cloud interactions by examining statistical 
relationships between aerosol properties and nearby low-altitude cloudiness using satellite data. 
The analysis of a global dataset of MODIS (Moderate Resolution Imaging Spectroradiometer) 
observations reveals that the positive correlation between cloudiness and aerosol optical depth 
(AOD) reported in earlier studies is strong throughout the globe and during both winter and 
summer. Typically, AOD is 30-50% higher on cloudier-than-average days than on less cloudy 
days. A combination of satellite observations and MERRA-2 global reanalysis data reveals that 
the correlation between cloud cover and AOD is strong for all aerosol types considered: sulfate, 


dust, carbon, and sea salt. 


The observations also indicate that in the presence of nearby clouds, aerosol size distributions 
tend to shift toward smaller particles over large regions of the Earth. This is consistent with a 
greater cloud-related increase in the AOD of fine mode than of coarse mode particles. The 
greater increase in fine mode AOD implies that the cloudiness-AOD correlation does not come 
predominantly from cloud detection uncertainties. Additionally, the results show that aerosol 
particle size increases near clouds even in regions where it decreases with increasing cloudiness. 
This suggests that the decrease with cloudiness comes mainly from changes in large-scale 
environment, rather than from clouds increasing the number or the size of fine mode aerosols. 
Finally, combining different aerosol retrieval algorithms demonstrated that quality assessment 
flags based on local variability can help identifying when the observed aerosol populations are 


affected by surrounding clouds. 
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1 Introduction 


Clouds, aerosols, and their interactions are among the largest sources of uncertainty in 
estimating human impacts on climate. As stated in the Intergovernmental Panel on Climate 
Change Sth Assessment Report (IPCC ARS, 2013) Chapter 7, “Clouds and aerosols continue to 
contribute the largest uncertainty to estimates and interpretations of the Earth’s changing 


energy budget ... the quantification ... of aerosol—cloud interactions continues to be a challenge.” 


Satellite observations have helped improve our understanding of aerosol-cloud 
interactions by providing observations on the relationships between aerosol and cloud properties. 
For example, analyzing a 4-year long global dataset of satellite observations, Loeb and Manalo- 
Smith (2005) found a positive correlation between cloud fraction (CF) and aerosol optical depth 
(AOD). Figure 1 shows that they found this relationship to be strong when using either NOAA 
or NASA satellite products. Other satellite studies also found positive correlations between CF 
and AOD [Ignatov et al., 2005; Kaufman et al., 2005; Zhang et al., 2005; Chand et al., 2012], or 
found that AOD increases in the vicinity of clouds [Loeb and Schuster, 2008; Tackett and Di 
Girolamo, 2009; Twohy et al., 2009; Varnai et al., 2013; Varnai and Marshak, 2014]. Similar 
tendencies were observed using ground-based and airborne measurements [Koren et al., 2007; Su 
et al., 2008; Ten Hoeve and Augustine, 2016] and also in modeling studies [Quaas et al., 2010; 
Grandey et al., 2013]. These studies identified several factors that likely contribute to the 
observed behaviors. For example, aerosol particles may swell in the humid air surrounding 
clouds [e.g., Twohy et al., 2009; Jeong and Li, 2010], or clouds may foster the generation of new 
or larger aerosol particles through chemical reactions or microphysical cloud processing 
[Kerkweg et al., 2003; Ervens et al., 2011, Eck et al., 2012]. On the other hand, higher aerosol 
contents may prolong the lifetime of clouds and thus increase cloud coverage in an area [e.g., 
Albrecht, 1989]. In some cases, remote sensing issues such as uncertainties in aerosol-cloud 
discrimination [Zhang et al., 2005] or the 3D radiative process of clouds scattering sunlight into 
nearby clear areas may also contribute [Marshak et al., 2008; Koren et al., 2009; Wen et al., 
2008; Varnai and Marshak, 2009]. Complementing the papers that focused on individual 
processes, other studies examined the relative importance of several factors in shaping the 
observed overall behaviors. For example, Jeong and Li [2010] provided insights on the way 
cloud contamination, aerosol swelling, and the formation of new aerosol particles affect the 


relationship between CF and AOD. 
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Figure 1. Relationship over global oceans between regional cloud fraction and mean aerosol 
optical depth. This figure is a reproduction of Fig. 6a in Loeb and Manalo-Smith [2005]. Both 
curves are based on data from the Terra satellite’s MODIS instrument, but they were obtained 
using two different cloud masking algorithms: a NOAA algorithm and the NASA algorithm used 


for producing the operational aerosol product. 


Despite the numerous insights provided by such studies, some important questions still 
remain. This paper examines two such questions: (7) How does the CF-AOD relationship vary 
with location and aerosol type? (ii) What is the relationship between cloudiness and aerosol 
particle size distributions? (We note that, in addition to being important for radiative transfer 
and other physical and chemical processes, particle size is also a key factor in the aerosol index 
that is becoming more commonly used in aerosol-cloud studies (e.g., Chen et al., 2014).) These 
questions are addressed through a statistical analysis of MODIS (Moderate Resolution Imaging 
Spectroradiometer) observations, sometimes in conjunction with MERRA-2 (Modern-Era 
Retrospective analysis for Research and Applications, Version 2) global reanalysis data on 
aerosol properties. First, Section 2 examines the relationship between cloud fraction and aerosol 
optical depth, then Section 3 examines cloud-related variations in aerosol particle size. Finally, 


Section 4 provides a brief summary. 
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2 Relationship between cloud fraction and aerosol optical depth 


Let us begin by analyzing the statistical relationship between CF and AOD (at 550 nm) 
using the MODIS Aqua Level 3 joint atmospheric product (Collection 6 version of MYD08 
product files, see Platnick et al. [2015]). This product provides daily average values of aerosol 
and cloud properties on a 1° by 1° horizontal grid covering the whole globe. (At low to mid 
latitudes this daily average comes from a single observation in the early afternoon.) Our study 
uses daytime average values for June-July-August (JJA) and December-January-February (DJF) 
from 2012 to 2014. We note that our analysis does not use observations that indicate complete 
cloud cover or completely clear sky in an area, due the lack of aerosol or cloud data on such 
days. Also, in order to focus on aerosol-cloud interactions for low clouds, data values were used 
only when the MODIS-estimated mean cloud top height of a partly cloudy grid box [e.g., Baum 
et al., 2012] was below 3 km. While this altitude limit eliminates the consideration of high 
clouds that are irrelevant to aerosols that occur mostly at lower altitudes, we note that some 
eliminated clouds that have high tops but low cloud base may have even greater influence on 


aerosols than the thinner, low-level clouds considered in this study. 


Figure 2a shows that the correlation between daily daytime average CF and AOD values 
for June, July, and August (JJA) is positive over much of the globe, with the exception of some 
desert areas. Panel b then shows that over most areas, AOD is substantially higher—typically by 
30-50%—on days when the cloud cover exceeds the local seasonal average than on days when 
cloud cover is below average. The results indicate that the global overall and regional 
relationships found in earlier studies truly extend over most of the globe (as opposed to coming 
either from a few dominating regions, or from some regions of the Earth being systematically 
both cloudier and more aerosol-laden than others). Also, the results are consistent with the 
results of Loeb and Schuster [2008], who found that cloud cover tended to be higher on more 
aerosol-laden days. We note that the results are similar for December, January, and February 
(DJF), though DJF data covers polar regions in the southern rather than northern hemisphere, due 


to seasonal changes in illumination and snow cover. 


In order to examine how the CF-AOD relationship depends on aerosol type, we combine 
the MODIS data with aerosol type information from the MERRA-2 global reanalysis. The 
reanalysis data is provided by the GEOS-5 (Goddard Earth Observing System Model, Version 5) 
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global model that incorporates the GOCART (Goddard Chemistry, Aerosol, Radiation, and 
Transport) aerosol model [Chin et al., 2002]. GOCART keeps track of five aerosol species: black 
and organic carbon, dust, sea salt, and sulfates. For simplicity, our analysis combines the black 
and organic carbon species into a single category. We note that while GOCART assimilates 
observations of aerosol amounts from MODIS, the distribution of various aerosol species comes 


from the model itself. 
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Figure 2. Geographical distribution of the relationship between cloud fraction and mean aerosol 
optical depth (at 550 nm) in JJA. (a) Map of CF-AOD correlation. (b) and (c) Maps of the 
relative and absolute differences between AOD values measured on days when cloud cover was 


higher or lower than the median cloud cover of each location. 
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Figure 3 shows the CF-AOD correlation derived from the MODIS data, but considering 
only those days when MERRA-2 indicated that a certain aerosol type dominated in the region 
(providing more than 50% of the total AOD in MERRA-2). Naturally, this means that 
observations from days with multiple evenly mixed aerosol types were not used. Moreover, the 
maps constructed separately for each aerosol type have large empty areas where the particular 
aerosol type never dominated. Overall, the results show a strong CF-AOD correlation for all 
aerosol types, even for the relatively less hygroscopic desert dust. This finding in consistent with 
the results in Gryspeerdt et al. [2016], which indicate that the AOD of various aerosol species 
calculated in the MACC global reanalysis increases systematically with CF. 


Correlation of CF & AOD 


Figure 3. Maps of local CF-AOD correlation considering for each location only those days in 
JJA when the MERRA-2 reanalysis indicated that a single aerosol type dominated in 1° by 1° 
area. (a) Days dominated by sulfate aerosols (b) Days dominated by carbon aerosols (c) Days 
dominated by desert dust (d) Days dominated by sea salt. Zero values (white) indicate locations 
where there were not enough days dominated by a particular aerosol type to calculate CF-AOD 


correlations. 
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To summarize this section, all observations either on a pixel-by-pixel basis (Level 2 data- 
product, e.g., Loeb and Manalo-Smith [2005]; Chand et al. [2012]; and Varnai and Marshak 
[2014]) or averaged over space and time (Level 3 product) provide a clear signature of a positive 
correlation between CF and AOD. The combination of satellite data with the global reanalysis 


model MERRA-? supports this correlation for all aerosol types considered. 


3 Cloud-related changes in aerosol particle size 


While all observations point to a similar statistical relationship between CF and AOD, 
this is not the case for the relationship between CF and particle size. To examine this 
relationship, Figure 4 shows global maps of the correlation between CF and aerosol Angstrom 
exponent (AE) in the MODIS Dark Target-Deep Blue combined product. In this product, AE 
data comes from the Deep Blue algorithm [Hsu et al., 2004] over land (for 412-470 nm), and 
from the Dark Target algorithm [Levy et al., 2013] over water (for 550-865 nm). Panel 4a 
reveals that in JJA, positive correlations dominate over much of the southern hemispheric 
oceans, but over northern hemispheric oceans there are large areas with both positive and 
negative correlations. Panel b shows that in DJF the two hemispheres are much more similar to 
each another, with positive but weaker correlations over most oceans. The figure also shows that 
large areas of both positive and negative correlations occur both over land and ocean. In other 
words, both the Dark Target and Deep Blue algorithms show that the effective aerosol particle 


size increases with cloudiness in some areas, and decreases in others. 


We note that reflectances are more sensitive to the size of fine mode particles at shorter 
wavelengths than at longer ones. This means that for aerosol populations dominated by fine 
mode, AEs involving short wavelengths will depend on the exact wavelengths used (e.g., Eck et 
al., 1999; Reid et al., 1999). Therefore, the wavelength difference between AEs produced over 
land and ocean may explain the jumps in CF-AE correlations that occur at many coastlines, for 


example at the East coasts of Asia and North America in Fig. 4. 


Let us mention that, in principle, the jump at coastlines can even help gain information on 


coastal aerosols if we assume that aerosol populations are similar on both sides of a coastline and 
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the accuracy of AE values [e.g., Sayer et al., 2013] is sufficient. This is because the shorter 
wavelengths used for Angstrom exponents over land are sensitive to the fine mode effective 
radius, while the longer wavelengths used for Angstrom exponents over ocean are sensitive to 
the fine/coarse mode aerosol fractions [Schuster et al., 2006]. Thus, having positive CF-AE 
correlations on land (suggesting smaller fine mode particles for higher CF) and negative 
correlations over ocean (suggesting lower coarse mode fractions for higher CF) could be a sign 
of increasing cloudiness being linked to decreasing values of fine mode effective radii and coarse 


mode fractions in coastal areas. 


We also note that while in Fig. 4, the correlations were calculated by giving all AE values 
equal weights, in some applications it is more appropriate to weigh AE values by the 
corresponding AOD. We found, however, that such weighting has little impact on the calculated 
correlations, and the maps based on AOD-weighted AE values (not shown) look almost identical 


to the maps in Fig. 4. 
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Figure 4. Map of the correlation between cloud fraction and Angstrom exponent in the MODIS 


dark target product over the ocean and the Deep Blue retrieval product over land. 


Seeking further insights into the CF-AE correlations, let us examine the CF-AOD 
correlations separately for the coarse mode and fine mode AOD values provided in the Dark 
Target product [Kleidman et al., 2005]. (The Deep Blue product does not provide separate AOD 
values for fine mode and coarse mode aerosols.) We note that in the MODIS Dark Target 
retrievals over ocean the assumed size of aerosol particles is fixed [Remer et al., 2005], and so 
variations in AE correspond to variations in the relative amount of coarse mode and fine mode 
aerosols. In good agreement with the results in Fig. 3 of Eck et al. [2010], simple Mie 
calculations show that a hypothetical increase of 0.1 in AE corresponds to a decrease of about 
0.05 in coarse mode fraction. This, in turn, implies a roughly 8% change in effective radius. To 
put this number in context, Mishchenko et al. [2004, 2007] argued that in order to reach an 
accuracy of 0.25 W/m’ in aerosol radiative forcing estimations, the uncertainty in aerosol particle 


size measurements should be less than 10%. 


A comparison of Figures 5a and 5b reveals that almost everywhere, CF is significantly 


more correlated with AOD for fine mode than for coarse mode aerosols. Figure 5c reveals that 
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over much of the open oceans, the AOD difference between cloudier and less cloudy days is 
larger for fine mode than for coarse mode. The main exceptions are regions with substantial dust 
transport from deserts and some areas near land. (The opposite trend in these areas may be due, 
in part, to a misidentification of some pixels with thick dust as cloudy, which may increase the 


estimated cloud cover in heavy-dust regions.) 


Let us point out that similar behaviors were also found in ground-based measurements 
[Eck et al., 2014]. For example, the analysis of AERONET sun photometer observations during 
the DISCOVER-AQ campaign over Maryland found that as clouds developed during the day, 
AOD increases were dominated by fine mode rather than for coarse mode aerosols (e.g., Figures 
5-7 in Eck et al. [2014]). The AERONET data also revealed that in some cases fine mode 
particle size increased with cloudiness as a result of humidification or cloud processing (their 
Fig. 6), while in other cases the fine mode shifted toward smaller sizes as a result of either new 
particle formation or humidification/cloud processing of very small particles (which made them 
large enough to become optically active) (their Fig. 7). We note that in all cases, the changes 
occurred at sizes much smaller than cloud droplets. This means that the basic findings cannot be 


attributed to cloud contamination, which should affect mostly coarse mode results. 


Similarly to this conclusion for AERONET data, a key implication of our results in Fig. 5 
is that over the large swaths of the oceans where AOD increases with CF more for fine mode 
than for coarse mode aerosols (i.e., red areas in Fig. 5c), cloud contamination is not the dominant 
reason for the CF-AOD correlations in MODIS data. We believe this based on the consideration 
that since undetected cloud drops are larger than most aerosols, they tend to have smaller AEs 
than bimodal aerosol populations. (The main exception would be areas where the fine mode is 
missing and coarse mode aerosols have slightly lower AEs than cloud droplets. However, this is 
not a typical case, as the average fine mode fraction is substantial in all oceanic regions [Remer 
et al., 2008] and, as it will be shown later in this paper, the mean aerosol AE values over most 
oceans are much higher than the near-zero or even slightly negative values we can expect for 
cloud droplets.) Given their small AEs, undetected cloud drops tend to reduce the AEs of aerosol 
populations. Because aerosol retrievals interpret a drop in AE as a sign of an increase in coarse 
mode fraction (e.g., Eck et al. [2010]), undetected cloud droplets tend to increase the retrieved 


coarse mode fraction. Therefore when cloud contamination is the dominant cause of CF-AOD 
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correlations, AOD must increase with CF more for coarse mode than fine mode—which is not 


the case in the red-colored areas in Fig. 5c. 


Naturally, this does not mean that cloud contamination does not occur (for example in 
partly cloudy pixels), or that its effect are not significant—only that in many areas it is not the 
dominant effect. The reason why it is not dominant is that data processing algorithms take extra 
care to minimize it. For example, the MODIS Dark Target algorithm shown in Figure 5 uses 
several steps to screen for clouds using various tests of brightness as well as spectral and spatial 
variability—and as an added precaution, it excludes the brightest 25% of pixels in an area even if 


they passed all tests [Martins et al., 2002; Remer et al., 2005]. 
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Figure 5. Map of relationships between CF and AOD separately for fine and coarse mode 
aerosols. Panels a and b show CF-AOD correlations for (a) Fine mode; (b) Coarse mode. Panel 
c shows where the AOD difference between cloudier and less cloudy days is larger (or smaller) 
for fine mode than for coarse mode aerosols. The yellow rectangle in Panel c highlights the area 


used for the analysis in Figures 6 & 7. All data is for JJA. 


While it is not yet certain why fine mode AOD varies more with CF than coarse mode 
AOD does, several factors may play a role, for example: (7) coarse mode aerosols—such as dust 
particles in the elevated Saharan Air Layer—may float well above boundary layer clouds and 
hence may not interact with them; (ii) fine mode aerosols are typically more hygroscopic and 
swell more in the humid air surrounding clouds; (iii) cloud processing or new particle formation 
may create small particles and/or more optically effective ones; (iv) sunlight scattered from 
clouds can bias satellite retrievals toward larger AOD values and smaller aerosol sizes through 
the 3D radiative process called bluing [Marshak et al., 2008]. To help better understand aerosol- 


cloud interactions, future studies are needed to evaluate the relative importance of such factors. 


One possibility for further analysis is to combine data from different sources. As a first 
step in this direction, we combine the MODIS Dark Target (DT) product with the MODIS Ocean 
Color product [Ahmad et al., 2010]. Namely, we examine the impact of clouds on AE by 
combining 3 km-resolution MODIS Dark Target AE values with 1 km-resolution Quality 
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Assessment (QA) flags that warn about possible cloud impacts in the MODIS Ocean Color 
product. To allow comparisons with earlier results, the analysis is based on the same MODIS 
Aqua observations as in Varnai and Marshak {2009, 2015]: September 14-29 in the 2002-2011 
period over the northeast Atlantic Ocean (45-50° N, 5-25°W). Let us point out, however, that 
based on Figs. 2, 4 and 5, this region appears to have relatively weak relationships between CF 
and AOD or AE. As a result, we expect that the tendencies discussed below are significantly 


stronger in many other regions of the Earth. 


Following the approach of Varnai and Marshak [2015], Figures 6 and 7 of this study 
show results for pixels that are 3 and 6 km away from the nearest low-level clouds (cloud tops 
below 3km). CF values in Figures 6 and 7 are calculated as the cloud cover of low-level clouds 
in a 41 by 41 km surrounding area. The figures show results for four different subsets of all 
available data. The black lines are based on all data where the Dark Target QA flag indicates 
high quality retrievals (DT QA > 1). The blue lines are for those pixels with DT QA > 1 where 
the Ocean Color QA flag “straylight” warns about potential impacts from clouds that occur 
within the surrounding 5 by 7 km _ region from a_ bright cloudy pixel 
(https://oceancolor.gsfc.nasa.gov/reprocessing/r2009/flags/). The red lines are for those pixels 
with DT QA > 1| where the ocean color QA flag “sstwarn” warns about the possibility of cloud 
contamination. (This flag is based largely on the variability of 250 m-resolution 0.86 um 
reflectances within 1 km-size areas and the variability of 1 km-resolution 11 um brightness 
temperatures in 3 km-size areas. For details, see https://oceancolor.gsfc.nasa.gov/atbd/sst/flag/.) 
Finally, the green lines are for those pixels with DT QA > 1 where neither the straylight nor the 
sstwarn flags give warnings. (Note that the sum of green, red and blue lines in Panel a can 
exceed the black line because some pixels have warnings by both the straylight and sstwarn 
flags. We note that the uncertainties due to annual variability (estimated from the spread of 


results for individual years) are around 0.0035 for AOD and 0.03 for AE values. 
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Figure 6. Impact of cloud-related Quality Assessment flags in the MODIS ocean color product 
on MODIS Dark Target aerosol data. (a) Number of pixels; (b) Aerosol Optical Depth (AOD); 
(c) Angstrom exponent. Figure is for September 14-29 in 2002-2011, for an area over the 
Northeast Atlantic Ocean (45-50°N, 5-25°W). Data is plotted for pixels that are 3 km away from 


the nearest low-level cloud. 
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Figure 7. Same as Figure 6, but for pixels that are 6 km away from the nearest low-level cloud. 


Figures 6a and 7a show that most pixels have cloud-related warnings 3 km away from 
clouds, but not 6 km away. (We note that most pixels with cloud contamination warning also 
have a near-cloud (that is, straylight) warning.) Figures 6b and 7b show that AOD does not 
change much with CF for pixels with no warnings about possible cloud effects, but it increases 
markedly for pixels with warnings. This suggests that the QA flags are largely successful in 
identifying cloud-affected areas (sometimes they might miss small thin clouds, or flag pixels that 
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are near not clouds but thick aerosol plumes). Figures 6c and 7c are more complicated to 
interpret. They show that fine mode aerosols dominated in the region, as all mean AE values are 
higher than the 0.75 value typically associated with 50% fine mode fraction of AOD at 500 nm 
[Eck et al., 2010]. Figures 6c and 7c also show that the (443-869 nm) AE provided in the 
MODIS ocean color product increases (that is, effective particle size decreases) with CF for 
pixels with no warnings. As suggested in Varnai and Marshak [2015], this may occur because 
small pollution particles may dominate aerosol populations when cloudy weather systems bring 
air from North America or Europe, whereas large dust or sea salt particles may dominate on less 
cloudy days when the air is coming from deserts or subtropical oceanic areas. In contrast, the 
effective particle size in pixels with cloud-related warnings does not change significantly with 
CF, as the decrease with CF seen for pixels with no warnings is offset by an increase related to 
cloud-related processes such as humidification. Particle size is markedly larger (AE is smaller) 
for pixels with cloud-related warnings than for pixels without warnings (especially 3 km away 
from clouds, see Fig. 6c), and this size difference increases steadily with the amount of clouds. 
This allows us to conclude that cloud-impacted (near-cloud) pixels and possibly cloud- 
contaminated pixels tend to contain larger particles (plus cloud droplets) than other nearby 
pixels. Finally, let us point out that the results in Figs. 6 and 7 imply that aerosol particle number 
concentration increases with CF both for pixels with and without cloud-related warnings 
(increasing AOD with steady particle size and steady AOD with decreasing particle size, 


respectively). 


Finally, Fig. 8 explores the relationship between the CF-AE correlation and AE itself. 
The map of mean AE values in Fig. 8a shows spatial patterns that are often the reverse of 
patterns for CF-AE correlations in Fig. 4a. This behavior is confirmed by Figure 8b, which 
shows that areas of small AE tend to have positive CF-AE correlations, whereas areas of large 
AE tend to have negative CF-AE correlations. This tendency can be explained by considering the 
hygroscopic swelling of fine mode particles, even if coarse mode particles remain unchanged. 
(For example, if coarse mode particles are less hygroscopic or float above the humid, partly 
cloudy boundary layer.) To demonstrate this, let us consider some simple Mie calculations for 
bimodal aerosol populations. The goal is to demonstrate that the same swelling of fine mode 
particles in broken cloud fields can have opposite effects on the overall AE, depending on the 


abundance of coarse mode particles. As an example, Figure 8c shows that if swelling increases 


16 


the median radius of fine mode particles from 0.07 to 0.12 um, the overall (550-869 nm) AE of a 
bimodal particle population will decrease along the dashed blue line in cases of large AE (low 
coarse mode fraction), but will increase along the solid red line in cases of small AE (high coarse 
mode fraction). We note that the AE of the bimodal population increases along the red line even 
though the AE of fine particles actually decreases. Specifically, as the fine mode median radius 
increases from 0.07 to 0.12 um and the coarse mode remains unchanged, the fine mode AE drops 
from 1.87 to 1.43 (not shown), the coarse mode AE remains constant at -0.13 (not shown), but 
the overall AE in Fig. 8c increases from 0.86 to 1.15. The overall AE increases because the 
swelling enhances the radiative impact of fine particles relative to the impact of coarse mode 
particles, and so the effective particle size of bimodal distributions shifts toward smaller sizes. 
(We note that even the blue curve would show an initial increase if swelling started at a smaller 
size: If fine particles started out so small that their light scattering was negligible, the overall AE 
would start out being the coarse mode AE. As fine particles then grew larger, the overall AE 
would increase to a value between the higher AE of fine-mode and lower AE of coarse mode.) 
The result that the same swelling (e.g., fine mode swelling from 0.07 to 0.12 um) can cause 
opposite changes in AE depending on the fine/coarse mode fraction, illustrates that the sign of 
AE variations by itself cannot reveal the direction of particle size changes. Consequently, the 
effect of clouds on aerosol particle size cannot always be adequately described by changes in 


Angstrom exponent alone. 


To summarize this section, the analysis of different MODIS aerosol products showed that 
over large regions of the Earth, effective aerosol particle size decreases with increasing 
cloudiness. This occurs because an increase in CF implies a larger AOD growth for fine mode 
than for coarse mode aerosols. The finding indicates that the observed correlations between CF 
and AOD do not come primarily from cloud detection uncertainties. Moreover, this section 
demonstrated a way in which cloud-related QA flags in the MODIS Ocean Color product can 
help identify MODIS aerosol data impacted by nearby clouds. Finally, the results revealed a 
duality in behaviors: Even when particle size does not increase with CF, it still increases in the 
proximity of clouds (e.g., Figs. 6 and 7 show larger particles at 3 km than at 6 km from clouds). 
This is possibly because variations in CF involve changes in the large-scale atmospheric 


environment (thus in aerosols), while the approach toward clouds occurs within the same large- 
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scale atmospheric environment and thus leads toward larger aerosol particles (as we see in both 


Ocean Color and Dark Target MODIS products). 
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Figure 8. Relationship between CF-AE correlation and particle size. (a) Map of Dark Target 
mean AE values for JJA; (b) Comparison of individual values in the maps in Figures 4a and 8a. 
The position of each dot along the x-axis represents the mean AE value of a grid point in Fig 8a, 
while the position along the y-axis represents the CF-AE correlation of the same grid point in 
Fig. 4a. (c) Change in the overall AE of bimodal aerosol distributions when the size of fine mode 
particles increases through hygroscopic swelling. The plot shows results from Mie calculations 
assuming Models #2 and #9 of the MODIS Dark Target (ocean) algorithm for the fine and coarse 
modes, respectively (Levy et al., 2013). As Table 2 of Remer et al. (2005) shows, the median 
radii are 0.06 um and 0.5 um for the “unswollen” fine and coarse modes, respectively, while the 
standard deviations of lognormal size distributions are 0.6 um and 0.8 um, and the effective radii 


are 0.15 and 2.5 um for the two modes, respectively. Hygroscopic swelling of fine mode 
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particles is simulated using the approach in Gasso et al. (2000). As fine particles swell from a 
median radius of 0.06 um to 0.28 um, their effective radius increases from 0.15 um to 0.7 um. 
The printed coarse mode fraction (CMF) values indicate what fraction of 550 nm optical 
thickness is due to coarse mode particles before any swelling. As fine mode particles swell, this 


fraction decreases. 


4 Summary 


The ultimate goal of this paper is to help better understand aerosol-cloud interactions—a 
leading cause of uncertainties in our estimates of human impacts on climate. To this end, it 
examines the statistical relationships between aerosol properties and the amount and distance of 


surrounding clouds in observations taken by the MODIS satellite instrument. 


First, the paper sought a more detailed picture of the positive correlation between cloud 
fraction and aerosol optical depth reported in earlier studies. Analyzing a global dataset of daily 
mean aerosol and cloud properties over a 1° by 1° latitude-longitude grid during three summers 
and winters, the study found that the positive correlation in earlier global statistics is strong 
throughout the globe in both winter and summer (as opposed to coming from a few dominating 
regions or from systematic differences between clouds and aerosols in different regions or 
seasons). Over much of the globe, aerosol optical depth (AOD) was found to be 30-50% higher 


on days that were cloudier than average than on days that were less cloudy than average. 


Combining MODIS observations with MERRA-2 global reanalysis data on aerosol type 
revealed that the correlation between cloud cover and AOD was strong for all considered aerosol 


types: sulfate, dust, carbon, and sea salt. 


After the initial focus on AOD, the paper examined how aerosol size distribution 
(characterized through the Angstrom Exponent) changes in the presence of nearby clouds. We 
found strong regional variations in the size distribution shifts that occur near clouds. Some large 
regions displayed shifts toward larger aerosols near clouds that is intuitively consistent with 
aerosol growth in the humid air surrounding clouds or the presence of undetected cloud particles, 
and is also consistent with CALIOP results in Yang et al. [2014]. Other large regions, however, 


clearly showed shifts toward smaller aerosols near clouds in three different operational MODIS 
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aerosol products: the Dark Target and Deep Blue products in this study, and in the Ocean Color 
product in Varnai and Marshak [2015]. The results also indicated that the quality assessment 
flags in the Ocean Color product can help us compare aerosol properties provided in the Dark 


Target product at pixels likely affected and by clouds and at pixels likely not affected. 


A more detailed look at the MODIS results revealed that when aerosol particle size 
distributions shift toward smaller particles in cloudier regions, this is due to a greater increase in 
the AOD of fine mode than coarse mode aerosols. The greater increase in fine mode AOD 
implies that the CF-AOD correlation discussed here and in earlier studies does not come 


predominantly from cloud contamination. 


Additionally, the results revealed a duality in behaviors: Even when aerosol particle size 
does not increase with CF, it still increases in the proximity of clouds. This is because variations 
in CF involve changes in the large-scale atmospheric environment and thus in aerosols, while the 
approach toward clouds occurs within the same large-scale atmospheric environment and leads 
toward larger aerosol particles. This implies that aerosol behaviors in partly cloudy regions are 
not dominated by a single factor and instead involve several factors: some change (increase or 
decrease) particle size throughout the cloudy regions, while others increase particle size in the 
vicinity of individual clouds. Such factors include hygroscopic swelling, meteorological 
covariation of cloudiness and aerosols, 3D radiative interactions (bluing), undetected cloud 
droplets, new particle formation through liquid phase chemical processes, and cloud processing 
that merges aerosol particles via the collision-coalescence and eventual evaporation of cloud 
droplets forming around them. The importance of each of these factors need to be evaluated in 


future studies. 


The results also revealed that the hygroscopic swelling of fine mode particles can either 
decrease or increase the overall Angstrom exponent of bimodal particle populations. The 
swelling of fine mode particles always decreases their own Angstrom exponent, but at the same 
time it can enhance the radiative impact of fine mode particles relative to coarse mode particles, 
if coarse mode aerosols are less hygroscopic or float above the humid, partly cloudy boundary 
layer. The greater relative importance of the fine mode can in turn lead to an increase in overall 


Angstrom exponent. This implies that the sign of Angstrom exponent variations does not 
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necessarily indicate the direction of particle size changes, and that the effect of clouds on aerosol 


particle size cannot always be adequately described by changes in Angstrom exponent alone. 


Although the methodology of data processing is also likely to impact the magnitude or in 
some cases even the sign of near-cloud particle size changes [e.g., Zgnatov et al., 2005], the 
finding of similar behaviors in three operational data products that use very different cloud 
detection, data selection, and aerosol retrieval algorithms suggests that the observed behaviors 
are not caused by data processing issues. We note, however, that CALIOP does not seem to 
show the diversity of behaviors observed by MODIS, and further study is needed to determine 
whether differences in cloud detection or data sampling may contribute to CALIOP consistently 


reporting larger particles near clouds [e.g., Yang et al., 2014]. 
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